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ABSTRACT 

In communication systems, efficient use of the spectram is an indis- 
pensable concern. Recently the use of compressed sensing for the 
purpose of estimating Orthogonal Frequency Division Multiplexing 
(OFDM) sparse multipath channels has been proposed to decrease 
the transmitted overhead in form of the pilot subcarriers which are 
essential for channel estimation. In this paper, we investigate the 
problem of deterministic pilot allocation in OFDM systems. The 
method is based on minimizing the coherence of the submatrix of 
the unitary Discrete Fourier Transform (DFT) matrix associated 
with the pilot subcarriers. Unlike the usual case of equidistant pilot 
subcarriers, we show that non-uniform patterns based on cyclic 
difference sets are optimal. In cases where there are no difference 
sets, we perform a greedy method for finding a suboptimal solution. 
We also investigate the performance of the recovery methods such 
as Orthogonal Matching Pursuit (OMP) and Iterative Method with 
Adaptive Thresholding (IMAT) for estimation of the channel taps. 

Index Terms: Sparse Channel, OFDM Channel Estimation, Pi- 
lot Allocation, Cyclic Difference Sets 

1. INTRODUCTION 

In wireless communications. Orthogonal Frequency Division Multi- 
plexing (OFDM) is a well-known solution for overcoming the prob- 
lem of multipath fading channels 1 1 1, |2|. However, this solution is 
effective only when the receiver is provided with tools to estimate 
the Channel Frequency Response (CFR). To this end, the transmitter 
should send some predefined data in a predefined order that the re- 
ceiver is a priori aware of. These predefined data are usually called 
pilots. 

There are two main approaches for inserting pilot data in 
OFDM signals. In block-type pilots, all the subcarriers in some 
OFDM blocks (the whole spectrum) are reserved as pilot tones. In 
comb-type pilot models, some predefined subcarriers in each block 
serve as pilots. Hence, CFR at these subcarriers can be estimated us- 
ing methods such as Least Square (LS) or Minimum Mean Square 
Error (MMSE). Now for estimating the CFR at non-pilot subcar- 
riers, interpolation methods ranging from simple linear or second 
order techniques 1 3 1 to time domain |4 1 and even more complex ap- 
proaches are used. It is clear that by decreasing the frequency gap 
between the adjacent pilot subcarriers, the performance of the inter- 
polation techniques improves. Therefore, the pilots are preferably 
put at equidistant subcarriers to provide uniformity. 

Considering the inherent sparsity in the impulse response of the 
wireless channels which is due to the sparse structure of the scatter- 
ing objects, it is possible to estimate the Channel Impulse Response 
(CIR) more accurately even from non-uniform pilot patterns. The 
common estimation techniques in this case are those introduced in 
the field of compressed sensing such as basis pursuit 1 5 1 and OMP 
171 . Unlike the interpolation case, equidistant pilot locations are not 
the best choices here. In [8], using the results of |5| for sparse sig- 
nal recovery, it is mentioned that uniformly random pilot location^] 
can provide the possibility of perfect channel reconstruction with 



overwhelming probability. Although this is an important theoreti- 
cal result, it is not practical. In this paper, we suggest a deterministic 
structure for the pilot locations in sparsity-based channel estimation 
methods which minimizes the inter-atom interference in DFT sub- 
matrices. Simulation results confirm the efficiency of the proposed 
pilot allocation method when greedy methods are used for chan- 
nel estimation. Also, we propose an iterative thresholding method 
for channel estimation which results in appropriate performance in 
time-variant frequency selective OFDM channels. 

2. PROBLEM STATEMENT 

In OFDM systems with comb-type pilot arrangement, ignoring the 
effects of Inter-Symbol Interference (ISI) and Inter-Carrier Interfer- 
ence (ICI), the received data at the fc''' subcarrier (1 <k<N)oi the 



OFDM frame can be formulated as: 



Y{n,k) = X{n,k) ■ H{n,k) + W {n,k), 



(1) 



where ^(n, k) is the transmitted OFDM symbol, H{n, k) is the chan- 
nel frequency response and W{n,k) is the AWGN noise.If !^ de- 
notes the set of all pilot indices, at a given pilot subcarriers kp !^ 
and using the least square method, the CFR can be estimated as: 
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(2) 



As explained earlier, conventional methods for estimation of 
the CFR at non-pilot subcarriers (given the noisy measurements at 
pilots) are interpolation-based techniques which require relatively 
high sampling rates (number of pilots) to produce acceptable Mean 
Squared Error (MSE). Also, the optimum structure of the pilot lo- 
cations for these techniques which minimizes the MSE of the es- 
timated channel, is the uniform distribution (equidistance) of the 
pilots in the spectrum. 

In the sparsity-based channel estimation methods, instead of 
finding the CFR, the goal is to estimate the inherently sparse CIR 
in each OFDM frame from limited number of noisy measurements 
of the CFR obtained at pilot locations. The estimated CIR is then, 
translated into the frequency domain by means of EFT which results 
in an estimation of the CFR that can be used for data equalization 
process. In these methods, we are dealing with the following system 
of equations: 



: Fp-h + Up, 



(3) 



where Fp is the DFT submatrix with Np 



' That means all possible choices of pilot indices ai'e equally likely 



rows associated 

with the pilot locations. Hp is the vector of LS-estimated CFR at 
pilot locations, h is the sparse CIR vector, and riy, is the vector of 
noise values. 

Generally, there are two main categories of sparsity-based 
methods to solve the set of equations presented in ([3]l. One ap- 
proach is to minimize the £] norm of h subject to (O, either directly 
or iteratively (such as SPGL |6|). Although the performance of 
such methods are considered among the bests, they are extremely 
slow for real-time implementation. The other approach which is 
considered in this paper, is to use fast greedy methods such as OMP 
which iteratively detect and estimate the location and value of the 



channel taps. These methods are usually faster than £i minimiza- 
tion techniques by orders of magnitude while they may fall short of 
performance. Our simulation results confirm that their performance 
is acceptable for the purpose of OFDM channel estimation. 

The main advantages of sparsity-based approaches can be cate- 
gorized into two parts: 

(1) Decreasing MSE: Generally, the purpose of using com- 
pressed sensing methods in solving a linear set of equations with 
the sparsity constraint is to achieve the Cramer Rao lower bound on 
MSE |9|. In extreme cases, the best estimator which achieves this 
lower bound is the structural LS estimator (usually called oracle es- 
timator) which exactly knows the location of nonzero taps (support) 
and estimates their corresponding values using LS 19|. However, in 
general, there is no information about the location of the nonzero 
coefficients of h at the receiver and the structural LS estimator is 
not realizable. Simulation results indicate that we can get close to 
this bound by using proper sparsity-based methods. 

(2) Reducing Overhead: Although the pilot subcarriers oc- 
cupy a fraction of the spectrum, they do not convey any data. By 
reducing the number of pilot subcarriers, we increase the utilization 
efficiency of the spectrum while we may degrade the performance 
of the channel estimation block. As mentioned in |8|, by consid- 
ering the sparsity of the CIR, it is possible to capture the necessary 
information in the frequency domain in fewer number of pilots. The 
results in |5| show that £i minimization technique almost perfectly 
reconstructs the sparse CIR from ([Sjl when the number of pilots is 
proportional to the number of channel taps. Furthermore, the re- 
construction performance is independent of the location and value 
of the taps; i.e., unlike the interpolation-based methods, the number 
of required pilot subcarriers does not depend on the delay spread 
and degree of frequency selectivity of the channel. 

3. ITERATIVE THRESHOLDING METHOD FOR SPARSE 
CHANNEL ESTIMATION 

In this section, we propose an Iterative Method with Adaptive 
Thresholding (IMAT)| 10| for the purpose of estimating the sparse 
CIR. In other words, we aim to identify non-zero channel taps and 
estimate their corresponding values using IMAT. 

In our general OFDM channel estimation problem presented in 
lO, our main goal is to estimate h from Hp given the fact that h 
has a few non-zero coefficients. To obtain an initial estimate, we 
multiply the sides of ([3} by Moore-Penrose pseudo-inverse of to 
find the solution with minimum i'2-norm: 

ho = F/ilp = Fp"''Fph + Fp^np. (4) 

Using the properties of the Moore-Penrose psuedo-inverse for un- 
derdetermind set of equations we have: 



:F/(FpF/)-' = iF/. 



(5) 



Now we can rewrite Q as: 



A' 



N 



(6) 



The elements of the N xN non-negative matrix G which is usually 
referred to as the distorting matrix, are given by: 



















G 




Threshold 






G 


a 













Threshold 



Figure 1: Block diagram of the IMAT method 



version of h; in general case where the columns are not orthogonal, 
ho is a distorted estimate of h. Now through a series of iterations 
and by employing the sparsity constraint, we try to improve this 
estimate. In each iteration, we perform one step of the iterative 
method studied in 1 11 1 followed by a thresholding operator: 



ht = A(ho-G-hj;_i) + h^_i, 
\k{i)\>l5e-"'< 

otherwise 



hd) = 



(8) 
(9) 



where A and k are the relaxation parameter and the iteration num- 
ber, respectively, G ■ h;(. is the output of the distorting operator to 
the input h^^ and ho is defined in l[6j. The steps {Sj are known to 
compensate for the non-orthogonality of the columns of Fy, while 
the thresholding operator takes the sparsity constraint into account. 
We employ adaptive thresholds in ^ which can be tuned through 
the parameters a and j3, and decrease exponentially with respect 
to the iteration number. The optimality of the exponential function 
in our method can be derived in a similar manner as in 1 13 |. The 
block diagram of the proposed channel estimation method is shown 
in Fig. [T] 

In the proposed method, the ideal but unrealistic case would be 
when there is no distortion; i.e., when G is a scaled version of the 
identity matrix. Although this never happens, it shows that when 
G is a good-enough approximation of the identity matrix, we can 
expect satisfying results. Thus, a given set of pilot locations is con- 
sidered as good if the off-diagonal elements of the matrix G are 
relatively small compared to the diagonal elements. In the next sec- 
tion, we will investigate the problem of selecting the DFT submatrix 
Fp which results in a suitable distorting matrix G. 

4. PILOT ALLOCATION BY MINIMIZING THE 
COHERENCE IN PARTIAL DFT MATRICES 

As mentioned before, the performance of a channel estimation 
block depends on both the reconstruction technique and the set of 
pilot locations. In this section, we will study the sub-optimum pi- 
lot locations when greedy sparsity-based methods such as OMP and 
the introduced IMAT are employed. 

4.1 Cyclic Difference Sets: Minimum Coherence 

To begin with, consider the following underdetermined set of equa- 
tions: 

ymxl = *mxnX„xl (10) 

where y is the observed vector, x is an 5-sparse vector (contains at 
most s non-zero elements) and $ is an m x n (in n) measurement 
matrix which in our case is the partial DFT matrix formed by 
selecting Np rows. In this section, we seek to find a proper location 
for pilots in each OFDM block, using the following definition. 



5'J ■ 




(7) 



where f,- represents the i'* column of F^. If the columns of F^ 
are orthogonal and there is no additive noise, ho will be a scaled 



Definition 1: The coherence of a measurement matrix €> e 
(^mxn jjjg fuaximum absolute cross-correlation between the nor- 
malized columns: 



\{(Pl ■ <Pk)\ 
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(11) 



Although the so called Restricted Isometry Property (RIP) f5] 
is the best known tool for characterizing the performance of a given 
matrix $ in sampling the sparse vectors, there is currently no poly- 
nomial time algorithm to check this property |12|. The common 
alternative for measuring how well a matrix preserves the informa- 
tion of the sparse vectors (x) in the produced samples (y) is the co- 
herence; the smaller the better. In addition, the performance of the 
greedy methods is more influenced by the coherence of the mea- 
surement matrix rather than its RIP order |7|. One of the well- 
known results demonstrates that the sparsity-based methods such as 
£i minimization and greedy methods, are guaranteed to perfectly 
recover the i'-sparse vectors when /i^ < ^ |7|. 

Returning to our main problem stated in ([3}, we aim to choose 
pilot indices in such a way that the coherence of the resulted mea- 
surement matrix, Fp, becomes as small as possible. Considering 
the unit-norm property of the elements of Fp, we have: 



max 

\<l.k<n 



l(f/,4)l 



A',, 



(12) 



If the pilot indices are = {Pi , . . .,P]\ip}, Fp becomes: 



it is obvious that 



, max {«,,}) > 



Np{N-Np) 
A'-I 



and the equality happens when 

g{i,{ad}) = --- = g{N-i,{aa}) 
which is valid only for 



Np{N-Np) 
N-l ' 



— — 
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(19) 



(20) 



(21) 



Hence, if there exists an index set ^ for which ai = • ■ • = ] 
happens, it is for sure the best possible choice for minimizing the 
coherence. Such index sets are already known as cyclic difference 
sets 1 14 1; unfortunately, the existence of difference sets are limited 
to some specific pairs {N,Np). 



^-jfP2(N-l) 



(13) 



According to the periodic structure of the DFT submatrix Fp, 
the inner product of f/ and f^, used in l ll2t only depend onr = k — l: 



max |(F,,F,+,)1 

l<r<A'-l 



max 

l<r<A'-l 



(14) 



1=1 



Here we aim to choose the set with = Np in order to 
minimize p.. For the simplicity of analysis, let us define /(x) = 

Y,iJli-^'- Hence, i fTil turns into: 



Mf„ 



max f(e 

l<r<N-\ '■ 



(15) 



Since we are interested in the modulus values of the function 
/(.) on the unit circle, instead of \f(x) \ it is simpler to work with 
= f{x)f*{x) = f(x).f{l/x). This shows that the optimum 
set 3^ which minimizes the coherence, is found by: 



opt 



= argmin max f{e '^)/(e-' '« 

6ja \<r<N-\ 



argmm max 

l<r<A'-l 



5^ ^e-^l^'-W-n-). (16) 



Z=lt=l 



If the set of cyclic differences of 3^ is defined a& S> = {Pi^ — 
Pl{modN) I 1 < l,k < Np; I ^ m} , and denotes the number of 
repetitions of the number 0<rf<A'— lin the set ^, we have: 



Y =Np+Y ade-jf'-" ^ g{r, {a,}). (17) 



Therefore, it is clear that we should look for the set of indices ^ 
which minimizes the maximum value of the function g(r, {a^/}) over 
all 1 < r < A'. Since 

N-l N-l 

L s{r,{ad})=Np{N-l)- ^ ad^Np(N-Np), (18) 

r=l d^l 



4.2 Greedy Coherence Minimization for Improper Pairs of A' 

and Np 

Cyclic difference sets described in section 14.11 are the optimal 
choices for the OFDM pilot locations with respect to the coherence 
criterion. In fact, if we make a DFT submatrix based on a cyclic 
difference set, the resultant matrix meets the Welch lower bound 
1 15|. In other words, not only is such a submatrix optimum among 
all DFT submatrices of its size, but also is the optimum code-book 
in the sense of minimum coherence among all the matrices with 
the same size. Nevertheless, for many pairs of N and Np, there 
is no cyclic difference set. Therefore we have to find proper in- 
dices for pilot allocation, using efficient search methods. In (16], it 



is suggested to use random exhaustive search among all 



\N„. 



DFT 



submatrices to find the one with the minimum coherence. As A' in- 
creases, the cardinality of the search space grows exponentially and 
the results of the random search in relatively small steps might not 
be satisfactory. Here we propose a greedy method to find suitable 
pilot index set. 

As stated in 14.11 it is important that the set of cyclic dif- 
ferences of the set of pilot indices ^ has equal number of 
repetitions (aj) for its different elements; i.e., the variance of the 
set of repetitions {0^}^ is equal to zero. In our greedy method, 
we choose Np pilot indices in the following Np stages: since 
rather than the exact value of the indices, their cyclic differ- 
ence are important, we initialize the index set by ^^'^ = {1} (1 
is arbitrary). The rest of the stages are summarized in the following: 

For the pilot index allocation: 

1 . Form all N — i+ I possible i-element subsets by adding an ele- 
ment to ^('-1): 



'Us 



5e{l,...,A'}\^('-i)} 



(22) 



2. For each ;-element set generated in step (I), form the set of 
cyclic differences and the set of repetitions ({aj}). 

3. Choose the set (or one of the possible sets) with the minimum 
variance in the elements of the respective repetition set ({a^})- 

4. If / < Np, go to step (I). 

5. SIMULATION RESULTS 

In order to give an insight toward theoretical results presented in 
sections [3] and |4l we conduct several computer simulations using 
MATLAB. The simulations are presented in two parts: 
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Figure 2: BER of different estimators at various SNRs 



5.1 IMAT Method for Sparse Channel Estimation Results 

The iterative method with adaptive thresholding presented in the 
block diagram of Fig. 1 is simulated in an OFDM system based 
on DVB-H standard with slight modifications. All the parameter 
specifications are presented in Table [T] For the channel, we have 
considered a Rayleight multipath fading channel with 4 significant 
nonzero taps at normalized (to carrier spacing) doppler frequency of 
1%; the average delay and power of the taps are presented in Table 

El 



Table 1 : simulation parameters 



Parameter 


Specifications 


Number of Subcarriers 


256 


Number of Pilots 


16(6.25%) 


OFDM Symbol Duration 


224{fls) 


Cyclic Prefix Length 


32(1/8) 



Table 2: Fading Channel Parameters 
Delay(^.s-) 1.7 3.5 5.2 11.3 
Power(dB) -2 -5 -7 



The IMAT method is compared with the linear interpolation 
method which estimates the channel at pilot frequencies using LS 
estimate (2) and then uses a linear interpolation function to estimate 
the CFR at data subcarriers. Also orthogonal matching pursuit is 
simulated as a proper sparse reconstruction method for channel es- 
timation. The obtained curves of the obtained Bit Error Rate (BER) 
and Symbol Error Rate (SER) shown in Fig. |2] and [3] indicate that 
the IMAT method outperforms the other competitors. 

5.2 Pilot Allocation in Sparsity-Based Estimation Methods 

In this part, we compare the MSE and perfect reconstruction per- 
centage in channel estimation for pilot allocation methods presented 
in this paper. For our simulations in this part, we generated a ran- 
dom 3-tap channel with varying fading parameters in each OFDM 
block and averaged the results over 5000 runs. Fig. |4] shows the 
MSE of the estimated channel for two different methods of pilot 
allocation. In the first scenario, the pilots are chosen uniformly at 
random for each block; in our proposed scheme, the pilots are ar- 
ranged according to a (73,9, 1) cyclic difference set and its cyclic 
shifts for different OFDM blocks. The Cramer-Rao lower bound is 
also presented in the figiu^es to give us a meaningful goal standard. 
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Figure 3: SER of different estimators at various SNRs 




SNR(dB) 

Figure 4: MSE of proposed and random pilot allocation methods 
for OMP reconstructuon. 



This bound is given by l9l : 

Ci?B = a2 trace (^(f„.a«Fp,a) , (23) 

where is the noise variance and Fp j^ is the submatrix of Fp ob- 
tained by keeping the columns corresponding to the channel taps. 
Cramer-Rao bound is a fair criterion to measure the quality of our 
pilot allocation method, since it is the MSE of an estimator that 
knows the exact location of the channel taps. Therefore, in OMP 
channel estimation, if we select pilot locations properly, the maxi- 
mum cross correlation between the columns of Fp becomes small. 
Hence, the columns of the resultant measurement matrix in (3) be- 
come less correlated which makes it easier to detect and estimate 
the CIR. Similarly for IMAT, decrease in the off-diagonal elements 
of the distorting matrix results in the decrease of the MSE values of 
the estimated channel (Fig. |5}. 

Finally, we compare successful channel recovery percentage in 
noiseless case for different pilot allocation methods. For this pur- 
pose, we have considered the OFDM communication with N = 256 
subcarriers and Np = 16 pilots. Since there is no cyclic difference 
set for this pair of an d Np, we have employed the proposed greedy 
search in section [431 to find a pattern for pilots with small coher- 
ence. The recovery percentage which is the percentage of exact 
channel recovery without error for various number of channel taps 
and OMP reconstruction method is presented in Fig. [6] 
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Figure 5: MSE of the proposed and random pilot allocation methods 
for IMAT reconstruction. 
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6. CONCLUSION 

In this paper, we investigated the problem of OFDM pilot allocation 
in sparsity-based channel estimation methods. First, we proposed 
an iterative method with adaptive thresholding (IMAT) which de- 
tects channel nonzero taps and their corresponding values iteratively 
for the purpose of OFDM channel estimation. As it was shown in 
the simulation results, this method outperforms typical interpola- 
tion methods such as Linear Interpolation and greedy algorithms in 
sparse channel estimation. We derived the optimum pilot location 
for greedy methods in sparse channel estimation, based on minimiz- 
ing the coherence in DFT submatrices. Simulation results show the 
improvement in the MSE of the estimated channel for our proposed 
pilot allocation method compared to uniformly random insertion of 
pilots. 



[5] E. Candes, J. Romberg, and T. Tao, "Stable signal recovery 
from incomplete and inaccurate measurements," , Comm. Pure 
Appl. Math., vol. 59,pp. 1207-1223, March 2006. 

[6] M. A. T. Figueiredo, R. D. Nowak, and S. J. Wright,, "Gradi- 
ent projection for sparse reconstruction: Application to com- 
pressed sensing and other inverse problems," , IEEE Journal on 
Selected Topics in Signal Proc, vol. 1, no. 4, pp. 586597, Dec. 
2007. 

[7] Tropp, J.A, "Greed is good: algorithmic results for sparse 
approximation," , Information Theory, IEEE Transactions on, 
vol.50, no. 10, pp. 2231- 2242, Oct. 2004. 

[8] Taubock, G.; Hlawatsch, F.; , "A compressed sensing technique 
for OFDM channel estimation in mobile environments: Exploit- 
ing channel sparsity for reducing pilots," Acoustics, Speech 
and Signal Processing, 2008. ICASSP 2008. IEEE International 
Conference on , vol., no., pp.2885-2888, March 31 2008-April 
4 2008. 

[9] Carbonelli, C; Vedantam, S.; Mitra, U.; , "Sparse Channel Es- 
timation with Zero Tap Detection," Wireless Communications, 
IEEE Transactions on , vol.6, no. 5, pp. 1743-1763, May 2007. 

[10] F. Marvasti, A. Amini, F. Hadadi, M. Soltanolkotabi, 
B. H. Khalaj, A. Aldroubi, S.Holm, S. Sanei, J. Cham- 
bers, "A Unified Approach to Sparse Signal Processing," , 
|http://arxiv e.org/abs/0902. 1 853v 1 , Feb 2009. 

[11] F. Marvasti, , "Nonuniform Sampling: Theory and Practice" , 
Kluwer Academic/Plenum Publishers, 2001. 

[12] Z. Ben-Haim, Y. C. Eldar, and M.Elad, "Coherence-based 
near oracle performance guarantees ," Wireless Communica- 
tions, IEEE Transactions on , vol.6, no. 5, pp. 1743-1763, May 
2007. 

[13] Soltanolkotabi, M.; Soltanahan, M.; Amini, A.; Marvasti, F; , 
"A practical sparse channel estimation for current OFDM stan- 
dards," Telecommunications, 2009. ICT '09. International Con- 
ference on , vol., no., pp.217-222, 25-27 May 2009. 

[14] Charles J. Colbourn and Jeffrey H. Dinitz; ,"CRC Handbook 
of Combinatorial Designs"CRC Press 1996. 

[15] Xia, P. and Zhou, S. and Giannakis, G.B.; ," Achieving the 
Welch bound with difference sets," Information Theory, IEEE 
Transactions on, vol.51, no.5, pp.1900-1907. May 2005. 

[16] Hochwald, B.M.; Marzetta, T.L.; Richardson, T.J.; Sweldens, 
W.; Urbanke, R.; ," Systematic design of unitary space-time 
constellations," Information Theory, IEEE Transactions on , 
vol.46, no.6, pp.1962-1973, Sep 2001. 



REFERENCES 

[1] ETSI EN 302 304, "Digital video broadcasting(dvb); transmis- 
sion system for handheld terminals (dvb-h)", November 2004. 

[2] ETSI TS 102 428 , "Digital audio broadcasting(dab);dmb 
video service; user application specification", June 2005. 

[3] M. Hsieh and C. Wei, "Channel estimation for OFDM systems 
based on comb-type pilot arrangement in frequency selective 
fading channels," IEEE Trans. Consumer Electron., vol. 44, 
no. 1, Feb. 1998. 

[4] R. Steele, "Mobile Radio Communications" . London, England: 
Pentech Press Limited, 1992. 



